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We analyze electron dynamics and topological properties of open double quantum dots (DQDs) 
driven by circularly polarized ac-magnetic fields. In particular we focus on the system symmetries 
which can be tuned by the ac-magnetic held. Remarkably, we show that in the electron spin 
resonance (ESR) configuration, where the magnetic fields in each dot oscillate with a phase difference 
of n, charge localization occurs giving rise to transport blocking at arbitrary intensities of the ac 
held. The conditions for charge localization are obtained by means of Floquet theory and related 
with quasienergies degeneracy. We also demonstrate that a topological phase transition can be 
induced in the adiabatic regime for a phase difference of n, either by tuning the coupling between 
dots or by modifying the intensity of the driving magnetic held. 



INTRODUCTION 

The study of topological features in driven systems is 
a fascinating topic, because of the emergence of non- 
trivial topologies once the system interacts with exter- 
nal driving fields^. Most of the topological features 
observed have been obtained under the assumption of 
adiabatic evolution, but also other interesting dynamical 
effects have been predicted in the non adiabatic regime 
0-0] • Among them, electron spin locking induced by 
linearly polarized ac magnetic fields [5| or spin blockade 
induced by ac magnetic fields [8] have been recently pro- 
posed. Therefore, investigation of different regimes in 
driven systems, and the interplay between ac-fields, spa- 
tial symmetries and, topological features in quantum sys- 
tems becomes a very promising field which could provide 
different mechanisms to manipulate electron charge and 
spin in confined systems. Spin qubits can be produced in 
confined systems as quantum dots, and their manipula- 
tion by driving the system with ac fields has been one of 
the most recently active topics, both experimentally [9- 
OJ] and theoretically 0, |, \±Mm ■ 

In this work we analyze the electron dynamics and 
current through a double quantum dot (DQD) coupled 
to contacts and driven by circularly polarized magnetic 
fields. We show that the driving field is able to induce 
charge localization (CL) and then transport blocking for 
arbitrary field intensity at electron spin resonance (ESR) 
configuration. This localization effect depends critically 
on the symmetries of the system and its appearance can 
also occur in larger size systems such as arrays of quan- 
tum dots, linear ions traps[l8||, optical lattices [l9|. and 
more generally, in tunnel coupled systems with a pseudo- 
spin degree of freedom. We will show that in our system, 
at resonance condition, CL occurs, if the ac magnetic 
field oscillates in phase opposition in the different dots. 
We discuss as well how CL is also reached off resonance 
for dots with different Zeeman energies at a particular 
value of the ac frequency, i.e. at the mean value of the 
Zeeman energy splittings of the dots. Coupling the sys- 
tem to leads and applying a dc voltage gives rise to elec- 




Figure 1: Schematic figure of a DQD coupled to circularly 
polarized magnetic fields oscillating in phase opposition 

tronic current through the DQD and then to characterize 
CL from the current. We analyze the driven current by 
means of the Density Matrix (DM) formalism within the 
framework of the Floquet theory. It allows to effectively 
deal with the secular terms, and to characterize the effect 
of decoherence in the current. 

In the last section we discuss how the topological prop- 
erties of the system are modified in the case where the 
phase of the ac magnetic field within each dot differs in 
7r. We demonstrate that a topological phase transition 
can be induced in the adiabatic regime by tuning the in- 
tcrdot coupling or by varying the intensity of the applied 
magnetic field. 

MODEL: 

Our model describes a DQD in a strong Coulomb 
blockade regime, with just up to one extra electron, cou- 
pled to leads. The states spanning our DQD Hilbert 
space are: {| T l>, | II), | tfl), I 4*)> |0)}. The DQD is 
coupled to a magnetic field B 1 = (B l x (t) , B l y (t) ,B\) by 
magnetic dipole interaction g/i^S • B (t) being i = L, R. 
The Zeeman splitting can be different for each dot due to 
different g-factors, differences in the nuclei polarization 
or inhomogeneities in the applied field. The ac-magnetic 
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field, that we assume circularly polarized, i.e. B x (t) 



Kc.x cos (ut 



i) and B l y (t) = B l ac y sin (ut + fc), have 



a parameter <^ characterizing the phase difference be- 
tween dots (FigQJ. 

The Hamiltonian reads: 



H(t) = H s {t) + H B + V 
H s (t) = H + H»{t) + Hl 



(1) 



being Hq = a e i^cr %d<r,i the dot Hamiltonian for the 
i dot, H*(t) = V „ /i; i .. ( ,(/.o,i.S,,. (fi = x,y) the 
coupling of the electronic spin with the external ac- 
field, Hf c = ^2iB\S Zy i the Zeeman splitting in the 
i dot, and H tLR = t LR J2a,i^ 3 'dljd aji the interdot 
tunneling Hamiltonian ( [lr = g = h = 1). The 
bath Hamiltonian H B = J2i,a,K e hK b i,K,<T b i,K,a repre- 
sents fermionic reservoirs with operators b\ K a /bi : K,cr 
for electrons with ck energy and K momentum, and 

a.K f V>i a.K^i^ + ^- c - ) cou pl es the DQD with 
the electron reservoirs. We consider Floquet theory be- 
cause the dynamics can be easily extracted from the 
quasienergy spectrum, where degeneracies can be related 
with localization. 

Symmetries: First we analyze the symmetries defined 
in the DQD, considering the parity symmetry (PS) n : 
{x — > —x} and the generalized parity symmetry (GPS) 
: {x —x, t — > t + T/2} , which usually play an im- 
portant role in driven systems, and classify the solutions 
according to a Z2 groupQ. Applying the parity trans- 
formation to Eq[TJ we obtain the condition for parity 
invariance: B z l = B z r, = e R = 0, <p = 0, and 
Bac a = Bac u ■ Note that if we consider the generalized 
parity operation, an extra minus sign coming from the 
time dependent term shows up, leading to a non invari- 
ant Hamiltonian. The way to obtain Ut invariance is 
by considering <j) = <f>2 — 4>i = ir, i.e. a tt difference of 
phase between the ac- fields in each dot. This leads to a 
GP invariant Hamiltonian, but the PS is broken. There- 
fore the difference of phase <\> switches between n and 
invariant systems. 

In the present paper we shall consider n and in- 
variant Hamiltonians, just by tuning <p to zero and tt 
respectively. We also fix B^ c x = B^ c y , and just consider 
the possibility of asymmetric Zeeman splittings (breaking 
both, PS and GPS). 

Also an internal symmetry in a single dot can be de- 
fined. In the linearly polarized discrete Z2 sym- 
metry is present due to the time dependence of the field 
intensity |_B(t)|. By contrary, in the circularly polar- 
ized case, the intensity \B\ is constant and the system 
presents a continuous symmetry for rotations along the 
z-axis Ur (6) = e~ l8Sz . Furthermore, in this last case a 
time translation i — >• t' is equivalent to a rotation along 
the z axis. This difference in the internal symmetry of 
the single dot in the presence of magnetic fields with dif- 



ferent polarizations leads to the lack of dynamical spin 
locking^ in the circularly polarized case. 

Master equation for the reduced density matrix: 
When a quantum system is coupled to a dissipative bath 
such as a fermionic reservoir, exchange of energy and in- 
formation occurs, leading to decoherence . In order to 
analyze how CL induced by circular ac-magnetic fields is 
affected by decoherence, and how this is reflected in the 
current, we develop a master equation within the frame- 
work of Floquet formalism, that allows to analyze the 
current through the system. The influence of phonons 
can be analyzed in a similar way, and an analogous effect 
in the current is expected. 

We start by considering the Liouville equation for the 
total DM in the interaction picture. After some manip- 
ulations, the master equation in the steady state for the 
reduced DM elements p a p becomes (see details in [2(1 21 1 
and Appendix^: 

(0)) 

(J, V 

" R-awV (°) P"P " R lp^H (°) ( 2 ) 

within the Markov approximation and Floquet basis, 
and being R\ ^ (0) the transition rates due to the cou- 
pling with the contacts, which couple the different Flo- 
quet states (see EqJ7]m Appendix^. In order to evalu- 
ate the current we calculate the time derivative of the 
number of electrons Ql = (— e) (N L + N dotiL ) in the 
left dot Ndot.L = J2a dl L do-,L, plus the left reservoir 
Nl = YIk (t"k l o°k,l,cj- Hence the average current is: 
(I) = (QL) = i([H(t),Q L }). 

RESULTS: 

We calculate the occupation probabilities in different 
states and the current for different field configurations. 
The results are compared with those of the closed system. 
We also consider different temperatures and analyze its 
effect on the decoherence. 

There are two cases to be considered, both in reso- 
nance condition B z j, = B z r = u. First, we deal with 
a n invariant Hamiltonian (cj) = 0). In this case, the 
quasienergies are not degenerate for all B ac (see FigJ^l 
left), but just isolated crossings appear for certain values 
of the field intensity. The time evolution of the occu- 
pation probabilities in the closed system (not shown in 
the paper) indicates that CL does not occur for this field 
configuration independently of the field intensity and fre- 
quency, and the electrons oscillate back and forth be- 
tween the dots. Our numerical calculation shows that 
all Floquet states become equally occupied in the steady 
state. At a finite voltage /xl — hr 3> T, B z , B ac , the cur- 
rent reaches a finite constant value (FigJ3J), which con- 
firms the lack of CL. 
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Figure 2: (Color online) Quasienergies vs B ac for = 
(left) and <j> — n (right), in resonance condition. The ab- 
sence (presence) of IIt symmetry for = (0 = 7r), leads to 
non-degenerate (doubly degenerate) quasienergies. The de- 
generacy, present for all B ac in the right figure, drives the 
system to CL. tLR = 1/10, and B z ,l = B z< r = lo = 1/2. 



If we consider a IIt symmetry invariant Hamiltonian 
(4> = 7r), the quasienergies become doubly degenerated 
for all B ac (Fig[21 right). The occupation probabilities 
show that the system is in CL regime (FigHJ. 
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Figure 3: Current vs t for the case of II symmetry (0 = 0) 
at resonance condition. B z l = B z r = lo = 0.5, B ac = 0.5, 
t LR = 0.1, T = 0.01, T = 10~ 3 and \i L - /i« > T,B z ,B ac . 
The inset shows the current vs time in the steady state {/) ~ 
2.5 x 10 -4 (3.5 x 10 _2 pA). All values are in energy units of 
lOOmT ( 5.788/ieV). (lo ~ 4.4GHz), the time unit scale is 
~ 0.23ns. 

The existence of IIt symmetry allows the Z2 classifica- 
tion of the quasienergies, and then the double degeneracy 
for all B ac , driving the system to CL regime. The cross- 
ings obtained between quasienergies with the same II or 
ILr symmetry are allowed by the Wigner-Von Neumann 
theorem because of the continuous symmetry of a sin- 
gle dot in a circularly polarized magnetic field, but they 
do not affect the occupation probabilities of the Floquet 
states. The symmetry difference between this case, and 
the case of a linearly polarized magnetic field (i.e. where 
the single dot symmetry is discrete) , is the reason for the 
absence of dynamical spin locking in the former. 

In the present case, i.e., for circularly driven fields we 
found a very interesting result: CL is obtained for all val- 
ues of B ac intensity. The reason is that the quasienergies 



manifold is degenerated for all B ac . Localization is im- 
proved by increasing the field amplitude, as we will show 
below. 
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Figure 4: Occupation probability in the left dot P vs time for 
Bac = 0.6, in presence of decoherence due to the leads. Note 
the spatial CL induced as a consequence of the quasienergy 
degeneracy between Floquet states with opposite IIt symme- 
try. The inset shows the time evolution (Eq[3} for the closed 
system. B Z>L = B Z , R = lo = 1/2, T = 0.01, T = 10~ 3 , 
tLR = 1/10 and = 7T, all in units of lOOmT (i.e. 5.788^ieV) 
for /xt - fiR » 0. 

It can be seen by direct comparison of FigQ] and the 
inset, that the effect of decoherence due to the leads, is to 
remove the coherent oscillations without breaking the lo- 
calization induced by the ac-field. In order to get a better 
insight of the effect of decoherence on CL we obtain ana- 
lytically the localization probability for the closed system 
(0 = ir and B Z<L = B zR = u): 

p(t) = |< Ti hiK*)>l 2 + IUzl*(*)>l a 



1 



(cos (t^f. 



At 2 

LR 



B 2 . 



<R 



(3) 



Eq(3]shows that CL depends on the ratio B ac jt^ R . From 
Fig H] we assume that, in presence of decoherence the lo- 
calization is given by the minimum value of the coherent 
oscillations in the closed system. We then can calcu- 
late the expected CL for the open system (Fig|5j) , which 
is given by: P m i„ = min (P(t)) = A 2 / (l + A 2 ), where 
A = B ac / (2t LR ). 

From the last results, we conclude that the existence 
of Ht symmetry drives the system to charge localization, 
as lo is in resonance with the Zeeman splitting. Once it 
is tuned off resonance, the quasienergies split, and local- 
ization is destroyed. 

Now we consider asymmetric Zeeman splittings. This 
is a very common situation in real experiments, where 
Overhauser fields or different g-factors make difficult to 
achieve a symmetric configuration. If we calculate the 
Floquet spectrum for the Hamiltonian H$ (t) ([1]), and 
look for the necessary condition for degeneracy between 
quasienergies, we obtain: lo = loav '■= (B z .l + B z .r) /2. 
Note that now IIt is always broken because of the asym- 
metric Zeeman splittings. For the frequency condition 
lo = loav, the quasienergies overlap as in Fig|2] (re-scaled 
to the average Zeeman splitting), leading to CL regime 
at = 7t. 
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Figure 5: P m in vs A in presence of coupling with the leads for 
the steady state (cf> = 7r). 

Finally Figj6] shows (I) /to for cj> = {0,7r} and asym- 
metric Zeeman splittings. 
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Figure 6: Current vs ui in the steady state for asymmetric 
Zeeman splittings. (I) drops to zero at u = ojav (i^Av = 0.75) 
for <j> = 7r, while for = the current is finite. _B z ,l = 0.5, 
B z ,r = l,B ac = 0.55, t L R = 0.1, and T = 0.01. Different 
temperatures have been considered, but their effect just affect 
to the transitory regime. 

This result demonstrates the critical difference in the 
current behavior for cp = {0, 7t} at luav, where in the last 
case,CL is induced and the current drops to zero. Now, 
we will consider the adiabatic regime and we show that 
also the phase difference is critical and determines the 
topology of the system. 

Topological phases in the adiabatic limit: 

In the adiabatic limit, spin systems coupled to mag- 
netic fields show a special feature. A non-trivial topolog- 
ical phase (characterized by integer numbers Z) arises for 
a single localized spin whenever the adiabaticity parame- 
ter v = u/2\B\ fulfills v < 1, being \B\ = ^ B 2 Z + B 2 ac ^. 
As the u) increases, the system is driven to the non- 
adiabatic regime, inducing a topological phase transition 
to a region where the topology is trivial. 

In the present setup (for simplicity we consider B z ^ — 
B z ,r), where the spin can oscillate between the two QDs, 
we characterize the new topological features due to the 
tunnel coupling. Considering both cases cp = {0,7r}, we 



calculate the first Chern number c\ that fully character- 
izes the topology of the system, and the adiabaticity pa- 
rameter it (see appendix for details). The result shows, 
in the = case, that the tunneling does not change the 
topological properties in the adiabatic limit (same result 
as for a localized spin) with c\ = ±1 for all values of tLR- 
In contrast, for the cj> = ir case, it is possible to induce 
a topological phase transition within the adiabatic limit 
just by tuning A = tLnJ\B\, This can be observed in the 
first Chern number for each state: 

c t (A) = ~ (1 + sign (1 - 2A)) (1, -1, 1, -1) 




(1,-1,1,-1) VA<l/2 
VA>l/2 



This result demonstrates that tLR — \B\/2 is a critical 
point for the phase transition, and that within the adia- 
batic regime we can induce a topological phase transition 
just by tuning the intensity of the field. 

In order to be in the adiabatic regime for the present 
configuration, it is required to be out of resonance oj <C 
B z . Therefore the previous analysis concerning topology 
for cf> — {0, 7r} applies for regimes where CL cannot oc- 
cur. It shows, as we commented previously, that different 
driving regimes in dynamical systems provide different 
ways to manipulate electrons. 

CONCLUSIONS: 

We analyze the electron charge dynamics and the tun- 
neling current through a DQD attached to contacts and 
driven by a circularly polarized magnetic field. We 
demonstrate both numerically and based in symmetry 
arguments that tuning the phase difference between the 
ac magnetic fields applied to the dots it is possible to 
achieve charge localization. This effect is robust and 
perdure in presence of decoherence, induced by the cou- 
pling with contacts. The effect of decoherence due to spin 
baths will be analyzed in a further work [5^. Finally, we 
demonstrate that out of resonance, in the adiabatic limit, 
we can drive our system to different topological phases 
by tuning the ratio ti,R/\B\. 

In summary, two level systems such as QDs coupled by 
tunneling and interacting with ac-magnetic fields show 
interesting features which depend drastically on the field 
polarization. The nature of the magnetic dipolar cou- 
pling —p, ■ B and the interplay between spatial and spin 
degrees of freedom coupled through the ac-magnetic field 
allow to achieve charge localization for arbitrary intensity 
of the ac magnetic field in the ESR regime. Also, in the 
adiabatic limit, we show that a topological phase tran- 
sition can be induced. The present results demonstrate 
the huge horizon of possibilities obtained by combining 
ac magnetic fields and QD arrays, with direct applica- 
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tion in quantum computation with spin qubits and also 
in topological quantum computation, although the lat- 
ter should be studied further for the case of non-abelian 
gauge theories. 
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2011-24331 and ITN, grant 234970 (EU) for financial sup- 
port. A. Gomez-Leon acknowledges JAE program. 

Master equation. Detailed calculation 

Following previous works 0, [2l|, we consider the Li- 
ouville equation for the density matrix and weak coupling 
with the leads. We calculate the reduced density matrix 
by tracing out the bath degrees of freedom. The assump- 
tions at this point are the approximate factorization of 
the density matrix between bath and system, and Markov 
approximation. The density matrix of the bath is given 
by a fermionic thermodynamical ensamble in equilibrium. 

We assume that our interaction Hamiltonian is linear 
in the system and lead operators: 

v = J2 = T E (/'L^- - K. :K d{.) ( 4 ) 

j K,i,a 

being T the coupling to the leads, FjK a lead opera- 
tor, and Qj a system operator (K is the momentum of 
the electron in the lead). We have considered multindex 
notation, being ±j — (if, j, <r, ±£), i = L,R, a =J [A, 
and £ = ± (this means annihilation/creation operator 
for the system). Identifying the coefficients we obtain 
Qi,a,+ = di j(T , Qi,o,- = —d>i t<T , Fx,i,o,+ = ^ Ki<y , and 

When we trace out the bath degrees of freedom, the 
terms involving the statistical ensambles can be written 
as: 

Gi (t - T ,t) = Tr s {FjFi (t - t, t) PB } (5) 

3,1 

G 2 (t-r,t)=J2 Tr B {Pi (t - r, t) F jPB } , 

3,1 

The calculation of these terms, also integrating in t, leads 
to: 

gi (muj - e^) + i = irDi (1 - n F (mu) - e^ v - m)) 

= g 2 (muj - e iiv )_ i 
g 2 {muj - = nDiUp (muj - - m) 

= gi (mu - e iiv )_ i 

being Di the density of states in the i lead, that we as- 
sumed to be constant in our regime, pi the chemical po- 
tential, and n F (E) = (l + e E ^ KBT ^)" the usual Fermi 
distribution function. 



The rates of the master equation re- 
quire the calculation of the matrix elements 
(<f> a (t) \QjQi (t - T,t) p s (t) \4>fs (t)), being \4> a {t)) 
the Floquet states of the closed system. Using Fourier 
series due to the periodicity of Floquet states, and 
defining Q {n) ]ap = ± f Q T dte*""*^* (t) (t)>. We 

obtain the full master equation in Floquet basis: 

(d t + is a0 )p(t) a0 = Y,( R U^(t) + Rl^At))p(t)^ (6) 

-K^ (*) p (t) vP - Rlp,v„ (*) p . 

where R 1 * (t) = £ M M „ (M) e"~ Mt , and the 
Fourier coefficients are: 

xgi {efjtv - mu) {)i 

The approximation considered in the paper 

Rl al,u. v ( M ) = R a/3,nu can be assumed if the 
quasienergies are far from the boundaries of the 
Brillouin zone. 



Adiabatic limit 

We define the adiabatic parameter following as: 
uj c : = sup { | Ann (t) | : n^m = 1, ...,4}, 

te [n,T 2 ] c [o,2tt/w] 

A mn = (m;t\ — \n;t) 

being ujq the level spacing to the first excited level 
and t the total time needed for a cycle. In the case 
(\> = — > wf =0 = 0/2 and the adiabaticity parameter is 
v^ =a = {l/2tLR (taking 0Jq~° = tLR)- The case <p = n is 
very different because the level splitting to the first ex- 
cited state is no longer given by the tunneling parameter. 
Hence, 

"t* = \^\B\ 2 +Ul R + At LR B z 

-^\B\* + ttl R -*t LR B z , 

U f=cj\B\/^\S\^+4^, 
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and 



= 7T _ 



(p—7T 



The adiabatic regime requires r < 1, and then small 
frequencies compared with all the other energy scales. 
For the adiabaticity condition we have assumed t^R < 
B z . 

In conclusion we state that for frequencies small 
enough (far from resonance), varying the tunneling pa- 
rameter, or the intensity of the magnetic field, a topolog- 
ical phase transition can be induced. 
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